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ID PERIODIC POTENTIALS 
WITH GAPS VANISHING AT k 



Abstract. Appearance of energy bands and gaps in the dispersion 
relations of a periodic potential is a standard feature of Quantum Me- 
chanics. We investigate the class of one-dimensional periodic poten- 
tials for which all gaps vanish at the center of the Brillouin zone. We 
characterize them through a necessary and sufficient condition. Po- 
tentials of the form we focus on arise in different fields of Physics, from 
supersymmetric Quantum Mechanics, to Korteweg-de Vries equation 
theory and classical diffusion problems. The O.D.E. counterpart to 
this problem is the characterisation of periodic potentials for which 
coexistence occur of linearly independent solutions of the correspond- 
ing Schrodinger equation (Hill's equation). This result is placed in 
perspective of the previous related results available in the literature. 
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1. Introduction 



A well-known achievement of Quantum Mechanics is the under- 
standing of the band structure of the energy spectrum for periodic 
potentials [1]. Since ever the characterisation of dispersion relations 
and gaps between bands of permitted energy has turned out to be 
of importance, due to their crucial role in the conductor-insulator 
properties of crystalline solids. 

In this scenario, the case of one-dimensional (ID) periodic po- 
tentials, despite its simplicity, is physically meaningful not only for 
pedagogical reasons, but also because it models real structures with 
a preferred direction (as nanotubes or nanowires). 

Our interest here is to study a sub-class of the very general problem 
of ^vanishing gaps\ focusing on real ID periodic potentials of the form 
W'^ + W' + vq, vq being a constant and W{x) changing sign after half 
a period. These potentials - and only these - turn out to have all 
the gaps vanishing at the centre of the Brillouin zone. 

Remarkably, one ends up with a (not necessarily periodic) ID po- 
tential of the form v{x) = W{x)'^ + W'{x) in several different fields of 
Physics, as in supersymmetric Quantum Mechanics [S], where spec- 
tral analysis combined with the formalism of SuSYQM has led to a 
large class of analytically solvable ID periodic potentials fI7 \ [50] , 
and peculiar features of supersymmetry breaking can be exploited 
|im [U]. Other examples are in the framework of the Korteweg-de 
Vries equation, via the Miura transform [14^ [5], and in the map- 
ping of a Fokker-Plank equation onto a quantum stationary problem 
[391 112] ) together with simulation methods related [3] . 

The outline of this paper is as follows. In Section [2] we state the 
physical problem and its mathematical formulation. Also, we men- 
tion some previous results on the criteria which make energy gaps 
disappear. In Section [3] we state and discuss our vanishing-gaps re- 
sults with some examples. Section 3] shows the proofs by means of 
elementary Quantum Mechanics and operator theory, as well as stan- 
dard O.D.E. and Complex Analysis theory. The appendix contains a 
more detailed review of the mathematics underlying the band struc- 
ture theory. 

2. Background Theory 

A Schrodinger-like particle in a periodic ID potential is described 
by the Hamiltonian 





4 



O. Zagordi and A. Michelangeli 



where v{x + 1) = v{x). A unit period for v, as well as units h = 
2m = 1, are assumed without loss of generality. We take v to be 
real and suitably regular, as specified later. This means that the first 
Brillouin zone is bounded by — vr, tt. 

According to the Floquet-Bloch theory, the eigenfunctions of H 
can be chosen of the form ipkix) = e^^^Ukix), with Uk{x + 1) = itfc(x), 
namely a plane wave e**^^ modulated by a periodic function of the 
same periodicity of v. The energy E{k) of the eigenfunction -0^, when 
plotted against k, gives the well-known band structure. Forbidden 
and permitted energies, bands filling, conductivity, and other related 
features can then be discussed in view of such dispersion relations 

Since we are interested in the vanishing of the gaps at the centre 
{k = 0) and at the edge {k = vr) of the Brillouin zone, we study 
the Hamiltonian (j2.1|) as the operator H{Q) or ^^(vr) defined on the 
domain 2?+ and P_ respectively, of the measurable functions on the 
interval [0, 1] that are square-summable together with their first two 
derivatives, and satisfy the boundary conditions 

V^(i) = V(o), i^'ii) = V^'W (2?+) , . 
V^(i) = -V(o), = -V^'W ^'-'^ 

One recovers the Bloch periodic/antiperiodic functions simply ex- 
tending any ij) G T)± by periodicity on M, and the familiar structure 
with bands and gaps as depicted in Fig. [1] with the customary nota- 
tions. Eigenfunctions at fc = have period equal to 1, while those at 
k = TT have period equal to 2. One band and its subsequent collapse 
into a unique band in the centre of the Brillouin zone at the energy 
E, iff £^ is a doubly degenerate eigenvalue of H{0). Similarly they 
collapse at the edge of the Brillouin zone iff is a doubly degenerate 
eigenvalue of H{'k). 

A more detailed description of this theory, which in the mathemat- 
ical literature is known as the Hill's equation theory, is reported in 
the appendix. 

Notice that the emergence of vanishing gaps is quite 'unusual', as 
for most potentials all gaps have a nonzero width [36^ I44j . In fact, 
in a sense that resembles the way that real numbers can be suitably 
approximated by rationals, any potential with vanishing gaps can 
be suitably approximated with a potential with all nonzero gaps. 
This establishes a sort of 'peculiarity condition' for potentials with 
vanishing gaps, which motivates the interest towards them. 

An example of ID periodic potential which is well known to exhibit 
some vanishing gaps in the dispersion relations is the square well 
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Figure 1. Typical dispersion relations for a ID periodic potential 



Kronig- Penney model [7] when certain conditions on the depth and 
length of the wells are met |331 I13j . Nevertheless, infinitely many 
nonzero gaps always occur. 

Through numerical analysis vanishing gaps have been found to 
occur for a hybrid of triangular wells separated by flat interstitial 
regions [33]. The emergence of vanishing gaps has been established 
by perturbative analysis on a smooth potential [38], dispelling the 
misleading idea that in some way only flat sections are associated 
with the vanishing phenomenon. 

The Whittaker-Hill potential, also known as the trigonometric Ra- 
zavy potential [Ml EH HOI IMl US] , is an example where all gaps except 
for a finite number of them, vanish only at the centre or only at the 
edge of the Brillouin zone. In the particular case n = we have the 
example given in equation (13.61) below. 

Only a finite number of nonzero gaps appear in the class of the 
Lame and associated Lame potentials [271 EOl EH HSl [Ml 12^: de- 
pending on the choice of the parameters entering their definition, 
they still exhibit a finite number of bound bands followed by an in- 
finite continuum band, and a well-studied pattern of vanishing gaps. 
Analogous results have been proved to hold for some complex-valued 
PT-invariant versions [501 [^[5T]. 

Beyond these examples, on the other hand, necessary conditions 
are known on a ID periodic potential v with a prescribed number of 
nonzero gaps in its dispersion relations. In particular [41j : 

(1) if no gaps are present, then f is a constant [H H?! [201 [21]; 
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(2) if precisely one gap occur, v is a Weierstrass elliptic function 
[211 El [23]; 

(3) if only finitely many gaps are present, then v is real analytic 
as a function on the reals [36l [TH [TT] ; 

(4) if v{x + a) = v{x) and all gaps at A; = vr are absent, then 
v{x + ^) =v{x) [i [221 [21]. 

In view of (4) above, we may assume that v has period 1 in the 
sense that 

1 = min{a > : v{x + a) = v{x)} 

(incidentally, this excludes v to be trivially a constant); if so, some 
gap at fc = vr must be open, whereas nothing is said at /c = 0. It is 
this question that we are facing in the following. 

3. Real periodic potentials vq + W{x)'^ + W'{x) 

WITH W{x + ^) = -W{x) 

We now come to the main object of our analysis. 
Theorem. Let v be a 1-periodic continuous potential. A necessary 
and sufficient condition for v to have all gaps vanish at the centre of 
the Brillouin zone is that 

v{x) =vo + W^{x) + W'{x) (3.1) 

for some constant vq and some differentiahle W changing sign after 
half a period: 

W{x + ]^) = -W{x). (3.2) 
Moreover, in terms of v, the function W is given by 



1 r+5 



W{x) = -- / v{0- / v{C)dC d^, (3.3) 

^ L Jo ^ 

thus it is determined only by the odd harmonic part of v: 

, v{x\ — v(x + i) , ^ 

W'{x) = ^ ^ , (3.4) 

whereas the constant vq is given by 

vo= [ [v{x) - W'^ix)] dx . (3.5) 
Jo 

Conditions 13. 1\) and ^3. 2\) uniquely fix W and hence vq to have the 
form and 113. 5|) respectively. Also, whenever ^3.1\) and 113. 2\) 

hold, then vq is the lowest energy in the dispersion relations of v, i.e., 
it is the ground state of the Hamiltonian H{0) with periodic boundary 
conditions. 
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In other words, the theorem states that given a ID periodic poten- 
tial V (with period 1) and W and vq as in (13. 3|) and (13. 5p respectively, 
then V — [vq + W"^ + W'] = if and only if all the gaps vanish at 
the centre of the Brillouin zone. In particular, when W is given as in 
()3.3p . VF^ + W' has a zero-energy ground state, and vq + W'^ + W' has 
ground state vq. Thus, if Egs is the ground state of any ID 1-periodic 

condition Egs — ^ Q necessarily implies that at least one gap 
exists at the centre of the Brillouin zone and the quantity Egs ~ 
can be seen as a 'measure' of such a non- vanishing phenomenon. 

We mention that a more detailed analysi^ shows that with a bit 
of standard (although non trivial) functional-analytic technicalities, 
the regularity of W can be considerably weakened. 

It is worth noticing that potentials characterised by the theorem 
above may or may not have some vanishing gaps at the edge of the 
Brillouin zone as well; nevertheless, as stated above, some of them 
must be necessarily nonzero, unless the potential is a constant. 

As an example let us concentrate on the very simple choice W{x) = 
sin27rx. The corresponding potential 

v{x) = W'^ix) + W'{x) = 2tt cos 2-Kx + sin^ 2-kx (3.6) 

has period 1, and is a particular case of the Razavy potential we men- 
tioned above. According to the theorem all gaps vanish at the centre 
of the Brillouin zone and the ground state is zero. The dispersion 
relations for v are plotted in Fig. [2l The emergence of nonzero gaps 
at the edge of the Brillouin zone is necessarily expected, our v being 
not a constant: indeed soniG opGn ^Q^ps cit k = vr are clearly visible. 
We underline that this very simple potential has only two Fourier 
components, making it interesting also in the field of optical lattices. 

When we slightly perturb W{x) = sin27rj; in such a way that the 
condition (13.20 is destroyed, e.g., by substituting 

W{x) 1 — > W{x) + er]{x) (3.7) 

where ry(x + ^) ^ — ri{x), the doubly degenerate eigenvalues at A; = 
split, and non-vanishing gaps appear separating the bands. Such a 
behaviour is reproduced in Fig. [3l 

As another example we point out the one-gap Lame potentials, 
which turn out to be a subclass of those we are dealing with, in that 
they show a single gap (the first one) at the edge of the Brillouin zone, 
whereas all gaps vanish at the centre. They have recently received 
new attention [H] since they optimise some key parameters in the 
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Figure 2. Centre: band structure for i;(a;) = 27r cos 27ra; + sin^ 27ra;. Left: 
close-up of the lowest energy band. Right: at A: = tt bands do not overlap. 
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Figure 3. Removal of the fc = vanishing gaps for the potential v 
corresponding to the perturbation W{x) = sin27rx i— > sin 27rx + e • sin 47ra; 



band structure and are of practical interest in the realization of quasi 
one-dimensional crystals. 

To summarise, we have seen that the entire class of potentials 
with all gaps vanishing at A; = is analytically characterised by the 
very simple formulas (|3.1|) and (|3.3|) . while the typical issue one can 
find in the literature is to identify specific potentials with a given 
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pattern of open and vanishing gaps (Kronig-Penney, Razavy, Lame, 
. . . ). Nevertheless, the formulas characterising this family might be 
used for band design purposes, or conversely to check whether a given 
potential belongs to this class. 

4. Proofs 

4.1. Proof of sufficiency. We want to make use of elementary op- 
erator theory, unlike the usual O.D.E. approach one can find in the 
literature (see appendix). For convenience let us rename the Hamil- 
tonian H{0), as defined in (|A.ip . as H := H{0). Once we set f o = 
the claim amounts to say that H is a positive operator, that its 
ground state is i^i =0, and that any other level -E > is doubly 
degenerate. Also, recall by (j2.2p and ()A.2p that the domain of H 
is characterised by the periodic boundary conditions 

^(1)=^(0), V'(1)=^'(0). (4.1) 

By means of the operator 



and its adjoint 



a := 4^ - W(x) (4.2) 
ax 



= --^ -W(x) (4.3) 
ax 



the Hamiltonian can be factorised as 

H = a+a. (4.4) 

Notice that W is regular enough to guarantee that ()4.4p makes sense 
without domain problems: both a and are only densely defined in 
L^[0, 1], nevertheless since W is differentiable, a maps into the 
domain of a\ so that a^a makes sense on the whole P+ (just as aa^ 
does). 

Of course factorisation (14. 4p means that H is a positive operator, 
hence all its eigenvalues are nonnegative. By direct inspection one can 
check that the ground state is i^i =0: indeed the (non normalised) 
state 

is annihilated by a, i.e., atpi = 0, since it solves the O.D.E. 

V''(x) - W{x)i;{x) = (4.6) 

and then also i/'i/'i = holds true. Ei is always nondegenerate, as 
we recall in the appendix. Notice, in addition, that the ground state 
wave function ipi can be chosen to be real. 



10 



O. Zagordi and A. Michelangeli 



We now want to exploit another factorisation, beyond the standard 
one (14. 4p . Let T be the operator of translation of half a period 

1. 



TiIj{x) = 'il){x + 



(4.7) 



defined with the natural periodicity in order that it becomes a unitary 
transformation on L^[0,1], i.e., 



1 



As a consequence of the assumption (|3.2p . the other operators 
transform under T as 

d _d_ 

dx ' 
-a , 



TWT^ 
TaT^ 



W, 

t 



-a 



THT^ = T{a)a)T^ = a a^ 
so that one can conveniently rewrite 

H = a)a = {Ta)\Ta) = {Ta){Ta)'^ 
and get the commutation relations 



[Ta,{Ta) 



[H,Ta] = H,{Ta 



(4.9) 
(4.10) 
(4.11) 

(4.12) 



It follows that any energy level is invariant under the action of Ta or 
(Ta)t. 

On P+ obviously T and T^^ square to unity; moreover Ta is anti- 



hermitian, i.e., {Ta) 
X £ [0, 1] 

(Ta^Tpix) = 



-Ta. Indeed, for any ip E and any 



a^T^ilj{x) 

-1p'{x-^)-W{x)^|Jix-^) 
^p'{x-^) + W{x)iP{x-^] 



(4.13) 



^'(x + i)-iy(x + i)V'(x + i) 

= —{Ta)'4j{x) . 

The anti-hermiticity of Ta implies that all its nonzero eigenvalues 
are pure imaginary numbers and that eigenfunctions belonging to 
distinct eigenvalues are orthogonal. 

Take now any -0 G !D_|_ and any A 7^ in iM such that Taip = Xip. 
Then Hip = \X\'^ip, because H = {Ta)^(Ta) = —(Ta)"^, that is, ip 
has energy E = |Ap. Taking all such possible choices of ^'s and 
A's, one recovers all the nonzero energy eigenvalues, because of the 
commutation relations (|4.12|) . 
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So, let Tmp = Xip: 

^'{^ + i) - W{x + ^)V'(x + ^) = AV(x) ; (4.14) 

by conjugation, since W{x) G M and A G iM, this is equivalent to 

i)'{x + ^) - + ^)^(x + ^) = AV^(x) = -AVi(x) , (4.15) 

V'(x) being the complex conjugate of V'(x). That is, TaV' = — AV'- 
Both ^ and have energy E = |Ap, but they are orthogonal, since 
A 7^ =^ A 7^ —A (distinct eigenvalues of To). Therefore, one has 
shown that the energy level E is doubly degenerate, and the conclu- 
sion must hold for any E > 0. □ 

Notice that conjugating the eigenvalue problem Ta^lJ = \ip turns 
out to be useful because it leads to the orthogonality of ip and V^. 
The same does not apply to the eigenvalue problem for H, in that 
Hil) = Etp do imply Hip = Eip, but one cannot argue that ip and ip 
are linearly independent. 

Concerning the possibility that some gaps (or all of them) vanish 
at the edge = tt of the Brillouin zone, it is clear that the scheme 
of the proof above does not apply. Indeed, at A; = vr the spectral 
analysis has to be performed on the Hamiltonian -ff(vr) now defined 
on the domain {antiperiodic boundary conditions); although the 
same factorisation ()4.1ip in terms of Ta still holds, when one tries to 
mimic (I4.13P one now gets 



(ra)t^(x) = atrt^(x) 

= -ip'{x - i) - W{x)'ip{x - i) 

= V'(x + ^)-t^(x + i)V(x + i) 

= {Ta)'p{x) 



(4.16) 



that is, Ta is hermitian on the antiperiodic functions. Then its 
eigenvalues are real numbers and conjugating Taip = Xip one gets 
Taip = Xip without being able to conclude whether ip^ip or not. 



4.2. Proof of necessity. Let us assume for the moment that v has 
mean zero: fmean := Jq vix) dx = 0. In the notation of the appendix, 

\F] IF] 

let ip[ and ■02 be the fundamental solutions and let D{E) be the 
discriminant of Hip = Eip, H being the Hamiltonian H{0) of the 
proof above, and of (|A.ip . and E being any complex number. 
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The object of crucial interest in this proof turns out to be the 
function 

f{E) ^ T^A2)^^^^ (4_17) 

^ ' l 2-D{E) ^ ' 

V E~Eq 

In fact the assumption of ah gaps vanishing at A; = translates into 

2- D{E) = M'^{E){E - Eo) (4.18) 

where all the zeroes of M{E) are simple and Eq is the lowest real 
root of D{E) = 2. Now, let Q^^\x) be the Green's function of the 
problem 

- -r^ + v{x) -E gl-^l (x) + 6{x) = (4.19) 

with periodic boundary conditions on [0,1]; by means of standard 
O.D.E. techniques, it is seen to be 

g[i.](,)^4^'(^)-V'fkx-l)^ ^4_20) 
^ ' 2-D{E) ^ ' 

Evaluating ^[^1 at mid-period (x = i) one gets 

^ ^2>~ M^E){E-Eo) ^ ' 

and since the boundary value problem ()4.19p is self-adjoint, such a 
Green's function can have only simple poles in the energy complex 
plane. Accordingly, by multiplication by Jv[[E)[E — Eq) 



G^-K-^)M{E)iE-E,,- ^^^^ 

.#k^)-#'(-^)__^(^) (4.22) 




D{E) 



-Eo 

must be an entire function. 

Our claim at this point is that / is also bounded (as E — > oo): 
then by Liouville's theorem it is identically constant. We enclose 
such claim in the following somehow technical Lemma. Its proof is 
postponed to the next subsection. Hereafter, the mainstream of the 
proof continues. 
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Denote by 



M{E) := 4'=l(i)-4'=I(-i) (4.23) 



the numerator and the denominator of / respectively. 
Lemma (uniform boundedness of / and estimates). The entire 
function f"^, and hence f, is uniformly bounded at infinity. Conse- 
quently it is identically constant, such a constant being 1 due to a 
direct evaluation. Further, in the region 

ne = {E eC : e ^ arg{E) s: 27r - e} (4.25) 

(e being any fixed small strictly positive number) the following asymp- 
totic estimates hold for the square of the numerator 



,,2. . 4sin2^ sin2^ .^/e^'^v^ 



c = 2 



E E^ \ ^5/2 

1 / /"^/^ \2 

vi^)+viO)-{l viOdC 



(4.26) 







and the square of the denominator 



4sin2^ sin2^ /e^-^^x , , 

M\E) = + AEo + O {-^) (4.27) 

where in both expansions each term is leading w.r.t. the subsequent, 
as E CO in Qe- 

Thanks to this Lemma the proof of necessity is completed as fol- 

lows. By comparison of the coefficients of the ( — ^^'^ j -terms in 
()4.26p and in (14.270 . one has c = AEq, whence 

2Eo =vi\) + viO) -1{JJ viO d^) ' . (4.28) 

Since the ground state Eq cannot change by any shift in x, this is the 
same as 

1 



2Eo = v{x + -)+vix)--(^J %(e)dc) 
= vix + \)+ vix) - ^ ( - 2Wix)^ ^ , 



(4.29) 
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where W is defined by ()3.3p - recall that we are now dealing with a 
zero-mean periodic v. Equivalently, 



v{x) + v{x + i) 



Eo + W'^{x). (4.30) 



2 

On the other side, definition (13. 3p clearly implies ()3.4p . namely 

ufx) — v(x + i) , , , , , 

^ ' ^-^ = W'ix), (4.31) 

so that altogether 

^(x) +t;(x + i) ^(x) -i;(x + i) 

= 2 + 2 (4.32) 

= Eo + W^{x) + VF'(x) . 

The constant Eq is recovered in terms of W by integrating over one 
period and taking into account that v{x) dx = W'{x) dx = 0: 

Eo = - [ W^{x)dx. (4.33) 
Jo 

So far V has had mean zero. For a generic v, the result above leads 
to V — V — fmean — -^'0 + + W , where W is given by the full form 
of ([33]) • Hence v = vq + W^ + W' with vq given by ([33]), and (f3l^ 
is proved. 

The uniqueness of is a consequence of the unique decomposition 
V = for any 1-periodic function v, where v±{x+'^) = ±v±{x). 

In fact, if -y = f + VF^ + W' = uq + U"^ + U' for some constants t^o, 
uq and some W, U changing sign after half a period, then necessarily 

= w = u' 

whence U = W + const, and by substitution into v the constant turns 
out to be zero. □ 

4.3. Proof of the Lemma. We proceed along the following steps. 

Step 1. f"^ is an entire function for which the following asymptotic 
estimate holds as E ^ oo: 

4sin2^ sin2^ ^e^^^- 
2 '4^0 — ^^ + 0( 



E ' " E^ \ ^5/2 
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with 



1/2 



(4.36) 



Notice that by no means this suffices to say that the ratio is 0(1): 
indeed the remainders are not meant to be necessarily subleading 
w.r.t. the sin^ and the E~'^ sin^ -^-terms, since the latter 
vanish in the sequences of points (2n7r)^. 

is entire because / is. The rest of step 1 is simply a consequence 
of plugging into ()4.17p the following asymptotic estimates available 
in the literature [Ml EI]: 



D{E) 



4sin^ 



smx 



E 



2 

E 



cos X 



^3/2 

"E 



E 



2E 



V(x) + 



(4.37) 
(4.38) 



+ 



smx 



[v{x)+v{0)-V^{x)] +0 



4^3/2 

i\/E G C and uniformly in x G [0, 1]), where 



£;2 



V{x) 



v{i) de . 



(4.39) 



These are extensively discussed in the appendix. We only remark that 
the continuity of v is all what is needed to truncate these expansions 
up to this order: to go higher to C(e'''"^^-E~"'/2)-terms one has to 
assume v to be sufficiently differentiable, which is not needed here. 
Also, we stress again that, until a given complex path — > oo is 
specified, it is not possible to identify the leading terms in (j4.37p and 
(08D . 

From (|4.38p on gets 



N{E) 



2sin^ 



+ 



sm 



E 



2E ^ ^2' 



v(4)) 



4^3/2 



(4.40) 



2 ^2' 



' g 2 2 
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Since v has period 1 and mean 0, then 



1/2 



v{0 



-1/2 



1-1/2 

/ v{Od^ = 0, 



V(- 



1, 



-1/2 



(4.41) 



v{i) de 

.1/2 

1/ 



/ V{i)di+ / ^;(e)dC = V(-), 
J-1/2 JO ^ 



whence 



2sin^ ^ sin^ 



E 



+ 



^2 



(4.42) 



Now, to square N{E) discarding the subleading terms recah that 



sm ■ 



^ e 



am 



cos ■ 



VSeC (4.43) 



(see, e.g., (jA.lSp in the appendix and the discussion thereafter) so 
that 



4sin2^ 
M\E) = ^ 

.1, 



+ c — z-r^ + O 



£;2 



^^(2)+^(0) 



with c = 2 
Analogously, from (|4.37p one gets 

M^{E) 



1/2 



2-D{E) _ ^^^^ ^ 



+ 



v{i) di 



^5/2 
2 ■ 



(4.44) 



4 sin 



2 -/E 



E 



+ 



Ell 



E 



(4.45) 
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SO that (I4.35P is achieved and the first step is completed. 

Step 2. Pick any e > small enough. Then is bounded in the 
closed region 0^ defined in ()4.25p . the bound depending on e: 

3Ce>0 : \f{E)\ ^ Ce yEen,. (4.46) 

In particular ()4.46p holds for every point on T^, the boundary of the 
angle C\0,£ with vertex at the origin, and 

lim = 1, E G^e. (4.47) 

E-*oo 

Here a shorter truncation in (j4.35p suffices, that is, 



4sin2^ sin2^ /e^"^v^\ 4sin2^ /e^'^^N 



E 



4sin2^ sin2^ ^e^"'^\ Asin^ ^ .g^mv^ 
Then 

4sin2^ ^e^mv^ 

hm f^(E)= lim ^ ' = 1 , (4.48) 



E 



^2 j 



because when £^ — > oo in the leading term both in the numerator 

and the denominator is E~^ svc? which dominates the O 

remainders. Indeed, when E £ fl^ and \E\ is large enough, (|4.43p 
actually becomes 



sin ^ e^-^ (4.49) 
(see, e.g., (jA.lSp in the appendix and the discussion thereafter), and 



VE J\E\ Be J\E\ e ^ . . 

V " 2^ T ^ 2^ sin - > , S E 0, . (4.50) 

To conclude this step, by continuity the entire function / can blow up 
neither at any finite point of nor at oo in 0,^ and ()4.46p is proved. 

Notice that if an arbitrarily small angle around the positive real 
axis was not cut off, then Um-^ would not necessarily increase up 

to +oo; furthermore the function sin^ has countably many zeroes 
in the real points En = (2n7r)2 and a priori it may not dominate the 

(d(^—f^^ remainder. 



Step 3. /2 has a finite growth order which does not exceed 1. 
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Recall ( |35j ) that an entire function /i : C — > C is said to have finite 
growth order p if 

In Mir) 

/9 := inf {/i > : M{r) < e''^} = limsup — < +00 (4.51) 

where M(r) := max^^^^j. \ h(z)\. Otherwise h is said to be of infinite 
order, and 

lim = +00 Vp > . (4.52) 

z—>oc gl^l 

We now see that our is not of infinite order. Ai{E) has only 
simple zeroes, that are all real; since / is analytic, these are also 
zeroes of M{E). From step 1 and V/i G M, 



e-l^l7'(i?) 



E 



is still an everywhere defined function, its numerator still vanishing 
whenever its denominator does; but now, due to (14.430 . as long as 
p > 1 



4sin2 ^ 



0' ^hwT;^^— (/^^l)- (4-53) 




This leads to 

namely ()4.52p fails to happen. That is, has finite growth order p. 
In view of (I4.5ip . we are sure that p is certainly dominated by the 
bound pb = 1- 

Step 4. is bounded also on C\Q,s; with the same bound as on O^. 
In fact the following facts turned out to hold: 

a) \P{E)\ Ce < +00 yE e Ts, the boundary of the angle 
C\Oe (step 2); 

b) is entire and its growth order does not exceed Pb = 1 (step 

3); 

c) the angle C\0,£ has arbitrarily small amplitude 2e, which can 

vr 

be taken less than — , namely vr. 
Ph 
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As a consequence of a standard corollary of the Phragmen-Lindelof 
theorem (see, e.g., theorem 9.12 in [35]), a)+b)+c) imply that 

\f{E)\^Ce yEeC\n,. (4.55) 



Conclusion of the proof. Fixed a small enough angle with ver- 
tex at the origin, centred around the positive real axis, and with 
amplitude 2e, the entire function turns out to be bounded both 
inside this angle, namely on C \ Og, and outside of it, namely on 
fig. That is, is bounded on the whole C-plane, and so must be / 
as well. Hence / is identically constant, by the Liouville's theorem. 
Such a constant can be evaluated, e.g., taking the limit £^ — > oo in 
some suitable way, as in (|4.48p . So = 1. □ 

The crucial role of analyticity of is remarkable. In fact, before 
applying Phragmen-Lindelof theorem for analytic functions, that is, 
without making use of analyticity, yet one has an apparently strik- 
ing control on it vanishes along any path ^ oo in C but an 
arbitrarily small angle centred at the origin around a given ray (the 
positive real axis). However this does not suffice to claim bounded- 
ness on the whole C. Neither it would suffice if in addition one knew 

to vanish at infinity also along the positive real axis, or along any 
ray emanating from the origin. A paradigmatic counterexample for 
this phenomenon is the smooth ^ M functions 

h{x,y) = {x^ +y^)e-^'=-y'^\ (4.56) 
Indeed for any e > small enough 
lim h{x,0) = 

x— >+oo 

N ^ V(x,?/)gR^ such that 

hm h(x,y) = n / \ r n i 

{x,y)-*oo 9 := a.ig{x,y) e [e,2TT - e\ 

because, when 6 G [e, 27r — e] (so that |sin0| ^ sine > 0) and 
r := Y^xM-y^ is large enough. 



(x - y2)2 = r^^sin^ Q _ cos ^ 



A ■ 4 

r sm e 



so 



^ h{x, y) < e-2'- ' > . 

r— >+oo 

Yet h diverges as on the curve x — y^ = 0. 



20 



O. Zagordi and A. Michelangeli 



Acknowledgements 

We are indebted to S. Baroni and G. Santoro for motivating the 
original interest to this topic. Helpful and critical discussions with 
R. Adami, V. Carnevale, G. Dell'Antonio, B. A. Dubrovin and G. Mor- 
chio are warmly acknowledged. 

Appendix. 

Spectral analysis of ID periodic potentials 

AND CONNECTIONS WITH HiLL'S EQUATION THEORY 

A.l — Spectral analysis. Degeneracy. As a starting point |4Hll9j 
one has to recognise that the space of wave functions on which the 
Hamiltonian (j2.1|) acts decomposes naturally into subspaces labelled 
by the boundary conditions at x = and x = 1. This way one has to 
consider the one-particle Hamiltonians 




(A.l) 



acting on the (dense) domain of self-adjointness Dk of the measurable 
functions on [0, 1] that are square-summable together with their first 
two derivatives, and satisfy the boundary conditions 

V'(l)=e*V(0) , V'(l) = eV(0). (A.2) 

(For convenience we rename P+ := and D_ := "Dn-) Thus the 
quantum-mechanical problem on M is rephrased in terms of the prob- 
lem of a Schrodinger-like particle on [0, 1] with boundary conditions 
labelled by k. 

The spectral analysis of H on the whole is the union of the spectral 
analyses of the -ff(/c)'s, and the following facts hold: 

• each H{k) has a purely discrete spectrum; 

• its eigenvalues Ei{k), E2{k), E^^k), ... are nondegenerate for 
any k G (0, vr); 

• -E'i(O) is nondegenerate as well; 

• /c En{k) is analytic in (0, vr) and continuous in [0, vr]; 

• for A; G [0, vr], A; I— > En{k) is monotone increasing for n odd 
and monotone decreasing for n even, and 

Si(0) < ^i(^) < ^2(vr) < ^2(0) < ^3(0) < ^3(vr) < • • • ; 

• H{k) and H{—h) are antiunitarily equivalent under ordinary 
complex conjugation, in particular their eigenvalues are iden- 
tical and their eigenfunctions are complex conjugates. 
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The dispersion relations k i— > En.{k) produce the famihar structure 
with bands and gaps (Fig. [T]). -En(") is the n-th band within the Bril- 
louin zone. Also, any ipk G T^k, when extended on M with boundary 
conditions hke (|A.2p on every interval [n, ?7- + 1], is the Bloch function 
tpk{x) = e'-^'^Ukix) with Uk{x + 1) = Uk{x). If 

_ J £'„(0) nodd o _j -E'n(vr) n odd 
"■~\^„(7r) neven' '^"■~\E„(0) n even ' ^ "^^ 

then [on, Pn] is the n-th band, On+i) is the n-th gap, i.e., the gap 
between n-th and the (n + l)-th band, and H has a purely absolutely 
continuous spectrum cr{H) = \J'^=i[an, Pn]- 

Eigenfunctions at the band edge k = have period equal to 1, 
satisfying the periodic boundary conditions. Eigenfunctions at the 
band edge k = n have period equal to 2, due to the antiperiodic 
boundary conditions. 

With this notation, the n-th gap vanishes iff Pn = Oin+i- this 
equivalently corresponds to En{0) = £'„_|_i(0) when n is even, and to 
En{T^) = -E'n+i(7r), when n is odd. Consequently one band and the 
following one collapse into a unique band at the centre of the Brillouin 
zone at the energy iff £J is a doubly degenerate eigenvalue of H{0). 
Similarly they collapse at the edge of the Brillouin zone iff is a 
doubly degenerate eigenvalue of H^n). 

A. 2 — Hill's equation theory. Main features. The dispersion 
relations k En{k) can be understood as branches of the parametri- 
sation [H [2] 

cos (Ve + 6(E)) , , 

of the variable E in terms of the parameter k, where d{E) and | t{E) \ 
are the phase and the modulus of the transmission amplitude t{E) 
of a free particle (i.e., a plane wave) of energy E incident on a single 
cell of the lattice described by v. Twice the left hand side of (1A.4|) 
is a quantity known as the discriminant D{E) of Hill's differential 
equation 

- ip" + V i/j = Eip (A.5) 

and this bridges the spectral analysis of ID periodic potentials to the 
Hill's equation theory [Ml EH EH 113 IS] • 
Such a discriminant, by definition, is 

D(£;):=VP(1) + (#1)'(1) (A.6) 
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Figure 4. Discriminant behaviour 



where E is now allowed to be complex, and V'l (x) and ■02 (^) ^'^'^ 
the so called fundamental solutions of the Hill's equation, that is, by 
definition, the solutions of the Cauchy problems 



4^1(0) 

(VfV(o) 



[E] 



[E] 



-(4 



4^1(0) 



1 







1 



(A.7) 



The discriminant turns out to be an entire function of the complex 
variable E with some crucial properties in dealing with the general so- 
lutions of (jA.Sp . Within this O.D.E. framework, what one commonly 
refers to as energy band, gap, band edge eigenvalues, vanishing gaps, 
translate respectively into the concepts of interval of stability, inter- 
val of instability, characteristic values, coexistence, with the features 
we sketch here below. 

• For real £"s, D{E) has a graph somewhat like that in Fig. HI 
In fact, the following asymptotic behaviour is known: 



£>(£;) = 2 cos + 0(^-3/2) E£R,E 
D{E) = 2 cosh ■ {l + o{l)) E£R,E 



+00 



-00 
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For complex E's, D{E) is an entire function of growth order 
^ and type 1 (|35j). 

• For any complex E, or real E such that > 2, then 
all nontrivial solutions of (jA.Sh are unbounded in (— oo, +cxo); 
unbounded solutions are called unstable. If E is real and 
\D{E)\ < 2 then all nontrivial solutions of ()A.5|) are bounded 
in (— oo,+cxd); bounded solutions are called stable. 

• D[E) = 2 has infinitely many real roots, which are single or at 
most double, increasing up to infinity, denoted in increasing 
order by 

ai,/32,a3,/34,a5,/36,--- — > +00 

and D{E) = —2 has infinitely many real roots, which are 
single or at most double, increasing up to infinity, denoted in 
increasing order by 

/3i,"2,/33,a4,/35,"6,--- — > +00 

so that whenever E belongs to any of the so-called stabil- 
ity intervals {an,Pn), then all the corresponding solutions of 
(1A.5|1 are stable, whereas if E belongs to any of the so-called 
instability intervals {[in-, On+i)) or even (—00, ai), then all the 
corresponding solutions of ()A.5|1 are unstable. 

• The n-th stability interval has a width (/?.„ — a„) < (2n— l)7r^ 
and — an) — (2n— l)7r^ — > as ?i ^ cxo. The asymptotic be- 
haviour of the ra-th instability interval width is (a„+i — — > 
and the positions of both edges On+ij l^n are asymptotically 
given by n^vr^ -|- v{x)dx . 

• All the solutions of (|A.5p satisfy the boundary conditions 

(|X2D 

^(x + 1) = e^^Xx) , + = e^*^V(^) 

or, equivalently, ii^ G R is an eigenvalue of H{k), if and only 
if D{E) = 2cosA;. In particular solutions V with D[E) = 2 
are 1-periodic, ^'s with D{E) = —2 are 2-periodic. This 
gives the correspondence between the bands/gaps and the 
stability/instability intervals. 

• Whenever D{E) = 2 has a double real root, say E = /?2n = 
a2n+i; then the 2n-th interval of instability (/?2n5 02n+i) dis- 
appears and there coexist two linearly independent 1-periodic 
solutions of (|A.5p with E = E (and therefore all the solutions 
are 1-periodic at that E). Analogously, coexistence happens 
whenever D{E) = —2 has a double real root E = P2n+i = 
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Oi2n+2'- there are two linearly independent 2-periodic solu- 
tions of ([Q]) with that E = E. 

Thus one sees that the vanishing gap phenomenon translates into 
the coexistence phenomenon for Hill's equation, associated to the 
double zeroes of 2 - D{E) = or 2 + D{E) = 0. 



A. 3 — Asymptotic estimates. We conclude this appendix dis- 
cussing the asymptotic expansions as — > oo of the fundamental so- 
lutions (|A.7p . Actually what we are doing in the following is to restate 
some classical but heterogeneous material [3^1 EH HZl [201 ED [Ml [E] 
in an organic unified perspective. 

The integral representation of the fundamental solutions is 



COS x^/E + 
sva.x\/E ^ 



E 



V-C/ 



(A.9) 



(these identities can be derived in a swift manner by Laplace trans- 
formation of (|A.7p ). whence, by iteration, 




sva.x^/E 



(A.IO) 



E 



+ 



E 

k=l 



cos x^/E 



smx 



+ 



COS X 



E 



EN 



+ 



smx 



E 



E^+"2 



Here Pj^^, QJ!' : [0,1] ^ R {i = 1,2; k = 1, . . . , N - 1) are the 

coefficients of '^"^^'^ and ^"^^/^ respectively in the expansion of the 

i-th fundamental solution. They are suitable continuous (and hence 
bounded) functions of x G [0,1], independent of the parameter E, 
and depending only on the choice of the potential v. Their existence 
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is guaranteed provided v is sufficiently differ entiable. In particular 
Pf\x) = 1. 

For instance, for real positive values of the parameter and if v 
is assumed to be only continuous (as in our main theorem), one finds 

r— sinx\/^ 

wi(x,E) = cosx\'E-\ ^^VuE 

2y/E 



^cos x\/^ 



Hx)-viO)-\v\x)]+o{^^ 
+ '-^^[v{x)+v{0)-\v\x)]+o{^) 



(A.ll) 



4^3/2 



with 



V{x) := /%(6de. (A.12) 
Jo 



The form of the remainders in ()A.11|) is just a shortcut to include 
two different types of remainders appearing in (|A.10p : both sinxV^ 
and cosx^/E do not exceed 1, since E G M4.. Incidentally, notice 
that plugging (|A.lip into the definition (|A.6p of the discriminant, 
immediately gives the first equation in (jA.Sp . 

The above expansions hold when ii^ — > 00 in the complex plane. 
The only difference with respect to the special real case is that in 
general the A^-th remainders can diverge (in modulus) as — > 00, as 
remarked in Section [4.31 

From fOO]) one deduces an analogous estimate for the quantity 
D{E), Indeed, plugging (IA.10|) into (jA.Gh and after some lengthy 
manipulations one gets 



E 

2-D{E) = 4sin2 



7V-1 



„;„ -/E „„„ v/E „■ 2 
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with 



0!k 

-27fc 



pW(i) + (p[2iy(i) 



Qi'ki) + (QL^')'(i) + n'H!i(i) 



[2] 



(A.14) 



As for the particular case of real £"s, it is customary to condense 
all the remainders of (jA.lOp and of (jA.lSp into a unique form. This 
is done by noticing that V a G C 



sm a 



cos a 



= ^ (cosham(2a) - cos9^e(2a)) ^ e'^'^^^a)! (A.15) 
|2 = i (cosham(2a) + cos9^e(2a)) ^ el^™(2a)| 



Thus, remainders can be written as ©(el 1^"") and ©(el ^mv^l^j-n) 

respectively. Also, the absolute value can be removed since x S [ 0, 1] 
and 3mVE ^ 0, due to the standard determination of the square root 
in the complex plane (if C 9 = \E\e^''<i'^^'> with ^ arg(£;) ^ 27r, 
then ^/E = y^e^^'s^^^/^ ^jt^ ^ arg(£;)/2 ^ vr, whence 3■m^/E ^ 
0). 

In the case of zero-mean potential, when suitably truncating (IA.10|) 
and ()A.13p one gets exactly the expansions ()4.37p and (I4.38|) used in 
our proof. As a further example, going to higher orders one would 
get 



2 - D{E) = 4sin2 



E sin ^ cos ^ 



2^3/2 
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+ 



/E 2E 



4^3/2 



1. 



v{x)+v{0) - -V{x) 



(A.17) 



+ 



^3 



[21 [21 

where the coefficients Q2 and P3 (according to the notation in 
()A.10p ) are given by 



1 ( V^jx) 
8l 3! 



v'^iOd^-Vix) [v{x)+v{0)] +v'{x) -v'{0)^ 



V(x) v\OdC + -[v\x) + v\0)] 
+ v{x)v{0) - v"{x) + v"{0) + V(x) [ f '(x) - v'{0) ' 
-^[v{x)+v{0)]. 



Notice that higher order derivatives of v occur. 
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